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Abstract 



Many known models, which generally use a factorization hypothesis, give 
a poor account of the decays B — > J/ip + ifW. Usually there is a free overall 
factor, which is fit to the data, so that tests of the models rely upon ratios. 
The models tend to give too much K* compared to K and too much transverse 
polarization compared to longitudinal. Our microscopic calculations, which 
use perturbative QCD, do well for both ratios. A microscopic calculation 
allows us to see how well factorization, heavy quark symmetry, and other 
features of various models are working. In the present case, agreement with 
the experimental ratios is dependent upon a breakdown of factorization for 
one of the amplitudes. 

PACS numbers: 13.20.He, 12.15.-y, 12.38.Bx 
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I. FACTORIZATION AND DATA 



Gourdin, Kamal, and Pham jl| and Aleksan et al. @ point out that many known mod- 
els all of which use the factorization hypothesis, give a poor account of the decays 
B -> J/i/j + K^. 

In most models there is an overall factor, generally called <22 H, which is fit to the data, 
so that tests of the models rely upon ratios of K* and K decays, and ratios of longitudinal 
and tranverse polarization in the J/ip + K* decays. The models tend to give too much K* 
compared to K and too much transverse polarization compared to longitudinal. 

Our microscopic calculations, which use perturbative QCD, do well for both ratios. Al- 
though the results have been published in some detail 0-0? the charmonium B decays 
deserve some further thought because of the present interest in them, and we will attempt 
to make self-contained at least the qualitative parts of our present remarks. A microscopic 
calculation allows us to see how well factorization, heavy quark symmetry, and other fea- 
tures of various models are working. In the present case, we find a serious breakdown of 
factorization for one of the amplitudes. 



More explicitly, the ratios under study and their experimental values are |10 



R = Br(B -+ Jji, + K) ~ 171 ± °' 40 ' (1) 
and, dividing the K* rate into a longitudinal polarization part T L and a transverse one T T , 

T L \ J 0.80 ±0.08 ±0.05 CLEO [10] 



Vj K * 1 0.66 ±0.10 ±0.10 CDF [11| ^ 



Our own results for the two ratios are 1.76 and 0.65, respectively (using Table IV of |7|]). 

Factorization implies that the decays depend upon a set of form factors for a current 
connecting B to K^*> . As a benchmark — yes, we know the K^*> is light — the relations that 



heavy quark symmetry [12-14] implies among the form factors lead to 



R = ml + f"* „ 2.38 (3) 



m 



B 



and 



r J k* m 2 R + 4m 2 



0.42, (4) 



J/1> 



which, although they do not agree with the data, are not bad as a representation of many 
of the models. For information, in each combination m\ + 4m jy^, the m\ comes from Tl 
and the 4m^ comes from IV. 

Why does our calculation work for the ratios when others do not? Most importantly, the 
factorization hypothesis fails. It does not fail uniformly. Its failure is significant only for the 
transverse polarization final state of B — > J/ip ± K*. In this amplitude the nonfactorizable 
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contributions are about half the size and opposite in sign to the factorizable ones, which has 
roughly the effect of turning the "4" into a "1" in the previous equations, and giving decent 
agreement with the (T LL /T) Kt data. 

Also, surprisingly in this context, we find the heavy quark symmetry symmetry predic- 
tions for the form factors of the factorizable parts of the amplitude work surprisingly well. 
One might expect significant differences due to a nonperturbative cause, namely that the 
wave functions or distribution amplitudes of the K and K* are different. A wave function 
difference at the origin is shown by data that gives unequal decay constants for the K and 
K*, and the shapes of the two wave functions are also different. We use the distribution 
amplitudes for K and K* worked out from QCD sum rules by Chernyak, Zhitnisky, and 
Zhitnitsky fl5 |. The upshot is that the form factors relative to the heavy quark symmetry 
predictions are good, and that small corrections and nonfactorizable contributions keep the 
two ratios from just being inverses of each other. 

Some details will now come. 



II. MORE DETAILED DISCUSSION 
Factorization 

We should state what factorization means in the context of B -»• J ftp + The 
relevant part of the effective Hamiltonian density is 

G 

H eff = -j=V cb V cs ary^c c^b, (5) 

where 7^ = 7 M (1 — 75) and the matrix element we want is, generically, 

M=(X^\H eff \B), (6) 

where ip = rj c , J/ip, ■ ■ ■■ 

The factorization hypothesis is that the charmed quarks which are created go into the 
ip and, except for the weak interaction vertices, are unconnected to other quarks in the 
process. We also assume that the outgoing charmed quarks in the ip have small transverse 
momentum relative to the direction of the ip. If the factorization hypothesis is valid, one 
can show 

M = -^-^V cb V cs (c 7m lcU (*7#>W, (7) 

that is, the matrix element is a product of two hadronic factors, 

(c7^lc)</, = (^|c7 mL c|0) (8) 

and 

(*y#W = (X\s^ L b\B). (9) 

Quantity N c is the number of colors, and to allow for effects of mixing with other operators 
one usually replaces (—1/N C ) with a constant 02 which is determined by data. 
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Non- factorization 

In Fig. [TJ, parts (a) and (d) correspond to the factorizable contributions, in the present 
context, and parts (b) and (c) to the nonfactorizable ones. 

Upon first view, it is easy to believe that the nonfactorizable contributions are small. The 
gluon couples to two oppositely colored quarks that are nearly at the same point because 
of the W -exchange. Indeed, the largest parts of diagrams (b) and (c) cancel each other and 
the subleading 0(qc), where qc is the gluon momentum, terms give the surviving result. 
The pieces of Figs. [I] (b) and (c) from one weak vertex, through the J / if) (whose polarization 
vector is £) including the gluon emmission vertex (7"), and to the other weak vertex, have 
numerators that sum to 

4mj/4l + 75) Ui v io - 4gi u (1 - 75), (10) 

which does go to zero for gluons of long wavelength, or qc going to zero. (The numerators 
of Figs, [l] (a) and (d) do not go to zero in the same limit.) However, the gluon momentum 
is not so small; in fact we argue that it is large enough that a perturbative calculation is 
plausibly valid. It supplies the momentum transfer needed by the light quark, which is of 
order A#, the part of the mass of the B meson carried by the light quark, which is about 
500 MeV or a few times Aq C d- 

However, for B — > J /if} + K and the longitudinal part of B — > J /if) + K*, there is further 
cancellation between the subleading parts of the two nonfactorizable diagrams. In contrast, 
they add for the transverse decay, so this nonfactorizable amplitude can get large. While 
the transverse B — > J /if) + K* does require chirality violations, the ensuing suppression is 
of 0(mj/^/mB), which is not a decisive factor. 



III. CLOSING REMARKS 

We have seen, in one explicit calculation, how nonfactorizable contributions to B — > 
J/ip + #M are significant and are crucial to giving agreement with data for the K* to K 
ratio and the tranverse to longitudinal polarization ratio. 

Ratios are used to test the models because in most models the overall rate is determined 
by a constant that is fit to the data. Our calculations using perturbative QCD also have 
trouble with the overall rate. At present, with the parameters we choose 0, we do a good 
job on non-color-suppressed decays such as B — > Dn but the overall rates are rather low 
compared to data for decays like B — > J /if) + K'*\ 

We believe the use of pQCD is valid for high recoil decays of the B. This has been given 



better support by Akhoury, Sterman, and Yao |16||who consider Sudakov effects in pQCD 
calculations of B decays and show that they suppress contributions from end point regions 
where use of pQCD would be questionable. 

It is possible that perturbative contributions are part but not all of what gives B 
J /if) + decay. In this case, the details of our remarks will only be part of something 
larger, but the scenario can well stand: the factorizable contributions to B — > J /if) + 
will not suffice to explain those decays, and nonfactorizable contributions will be crucial. 
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APPENDIX A: FACTORIZATION HYPOTHESIS 



Here is a perturbative proof that the factorization hypothesis leads to the factored form 
of the amplitude. Take the J /if) as an example. Factorization allows us to isolate the if) 
piece, 



sa^{l - l5 )(J/iP\ Cl c,\0)Y(l - 75)6,- 

fj/i> $ij £{4 + m J/i>) 



SiJ^l - 7s) 

1 



2y/W c y/N c y/2 



y(l-7fe)6i 



(Al) 



where q and £ are the momentum and polarization vectors of the J/ if) and i and j are color 
indices. One needs to recognize 



(J /ip\ey a c\0) = V^fj/ipmj/^a 



(A2) 



to complete the proof. 



APPENDIX B: FACTORIZABLE AMPLITUDES 

It is of some interest to see the behavior of the factorizable parts of the amplitudes. In 
practice, diagram 0(d) is quite small and we will give the results from Fig. |l](a). We use 



form factors f±, a±, g, and / as defined in ref. ||14|| . Neither /_ nor a_ enters when the decay 
involves the J /if). If desired, one may convert to ref. |§ definitions by (q 2 <-> here), 

Fi = U 

q 2 

^0 ~~ — 2 rf- + f+ 

m B — m K 
V = -(m B + m K *)g 
M = {m B + m^y 1 f 
A 2 = —(m B + m K *)a + 

A = {2m K ,)- 1 {f + {ml - m 2 K ,)a + + q 2 a_) (Bl) 

We keep the explicit m K (*) mass terms although they turn out to have small effect (unless 
necessary in a definition). 

Each form factor can be written like 

U = Bfi Vl ^ l -^ a + ^\ (B2) 
Jo yi — r — ii] 



where 



3e B (m 2 B + m 2 RM - m 2 ^) 2 



for e B = A B /m B and 
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Quantity 0* is the distribution amplitude of the kaon with the asyptotic form factored out 
(so that 4>k = 1 if we wish to use the asymptotic form). For form factor / + , 

a = a(/+) = m B (m B + m K ) - e B m B (m B - m K ) 
b = Kf+) = m \~ 2m B m K - rrij /f . (B5) 
For the other form factors, 

a{g) = m B (l + e B ) 

b(g) = -m K * (B6) 

and 

a(f) = 2m B m K *(l - 2e B ) 

+{m 2 B + m 2 K * - m 2 jj^)m B {l + e B ) 
b(f) = m K * (m| + m 2 K , - 4m B m K * - m 2 ^). (B7) 

It happens that the heavy quark symmetry results a± = ±g are obeyed exactly 
The integrals can be done analytically, using 

AT N-l , 

for N > 1 and 

/o = in + In ■ (B9) 

For asymptotic wave functions 



\9\ 



m B + m K ' m ' (m B - m K *) 2 - m 2 ^ 

= 0.99/jf : 1.03/jf. : 0.99/jc.. (BIO) 

Heavy quark symmetry predicts 1:1:1. Of course, strict heavy quark symmetry also predicts 
fx*/ Ik = 1 whereas the experimental result is 1.67/ 7r /1.22/ 7r w 1.37. 

One should perhaps not use the asymptotic distribution amplitudes for the kaons. A 
common form for representing the distribution amplitude is 

0( yi ) = 5/3(2 yi -l) 2 + (l-/3), (Bll) 

where y\ is the momentum fraction carried by the nonstrange quark and (3 is the fraction of 
the distribution amplitude that is not asymptotic. The QCD sum rule results of Chernyak, 
Zhitnitsky, and Zhitnitsky, lead to (3 = 0.6 for the K and 0.1 for the K*. One then gets 

1.38/* : 1.03/** : 1.03/*. (B12) 

for the same ratio as eqn. ( |B10| ) (with the same overall constants). This is stunningly close 
to the heavy quark symmetry result! 
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FIGURES 



FIG. 1. Lowest order perturbation theory diagrams for B decays involving charmonium. (a) 
and (d) are factorizable, (b) and (c) are nonfactorizable. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9409261vl 



